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^\ • Abstract 



We describe the generalization of spherical field theory to other modal 
expansion methods. The main approach remains the same, to reduce a 
j^^ . d-dimensional field theory into a set of coupled one-dimensional systems. 

^\ ' The method we discuss here uses an expansion with respect to periodic- 

O^ i box modes. We apply the method to </>* theory in two dimensions and 

compute the critical coupling and critical exponents. We compare with 
lattice results and predictions via universality and the two-dimensional 
fsing model. [PACS numbers: 05.f0-a, 05.50-fq, ff.fO.-z, ff.30.Qc, 
f2.38Lg] 



1 Introduction 



Recently spherical field theory has been introduced as a non-perturbative method 
C^ ' for studying quantum field theory [Q . The starting point of this approach is to 

decompose field configurations in a d-dimensional Euchdean functional integral 
as linear combinations of spherical partial waves. Regarding each partial wave 
as a distinct field in a new one-dimensional system, the functional integral is 
rewritten as a time-evolution equation, with radial distance serving as the pa- 
rameter of time. The core idea of spherical field theory is to reduce quantum 
field theory to a set of coupled quantum-mechanical systems. The technique 
used is partial wave decomposition, but this can easily be generalized to other 
modal expansions, fnstead of concentric spheres and partial waves, we might 
instead consider a generic smooth one-parameter fanuly of (d — l)-dimensional 
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manifolds (disjoint and compact) and basis functions defined over each mani- 
fold. There are clearly many possibilities and we can optimize our expansion 
scheme to suit the specific problem at hand. 

One important and convenient feature of spherical field theory is that it 
eliminates the need to compactify space — the spherical quantization surface 
is already compact. This is useful for studying phenomena which might be in- 
fluenced by our choice of boundary conditions, such as topological excitations. 
Another interesting example arises in the process of quantization on noncom- 
mutative geometries. Spherical field methods have recently been adapted to 
study quantum field theory on the noncommutative plane |g]. In many in- 
stances, however, maintaining exact translational invariance is of greater value 
than non-compactness or exact rotational invariance. In that case one could 
consider a system whose spatial dimensions have been compactified to form a 
periodic box. The next step would be to expand field configurations using 
free modes of the box and evolve in time. We will refer to this arrangement as 
periodic-mode field theory or, more simply, periodic field theory.R 

In periodic field theory the Hamiltonian is time independent and linear mo- 
mentum is exactly conserved. In this work we describe the basic features of 
periodic field theory and use it to analyze spontaneous symmetry breaking in 
Euclidean two-dimensional 0"* theory. We use the method of diffusion Monte 
Carlo to simulate the dynamics of the theory. The techniques discussed here 
have several advantages over conventional Euclidean lattice Monte Carlo meth- 
ods. One is that periodic boundary conditions are not imposed on the time 
variable, making it easier to determine the particle mass from an exponential 
decay fit. Another is that the zero- momentum mode is regarded as a single 
degree of freedom (rather than a collective mode on the lattice), which provides 
a simpler description of vacuum expectation values and symmetry breaking. 
Other advantages arise in systems not considered here, such as the absence of 
fermion doubling and extensions to Minkowski space using non-stochastic com- 
putational methods. 

The organization of this paper is as follows. We begin with a derivation 
of periodic field theory. We then analyze two-dimensional cj)'^ theory and its 
corresponding periodic-field Hamiltonian using diffusion Monte Carlo methods. 
We compute the critical coupling at which — > — refiection symmetry is 
broken, and determine the critical exponents f and /3.Q We find that our value 
of the critical coupling is in agreement with recent lattice results , and our 
values for the critical exponents are consistent with predictions via universality 
and the two-dimensional Ising model. 



^ Periodic field theory could be viewed as a hybrid of the Hamiltonian and momentum-space 
lattice formalisms. This combination, however, has not been discussed in the literature. 

^We are using standard notation, u is associated with the inverse correlation length, and 
P corresponds with the behavior of the vacuum expectation value of (f>. 



2 Periodic fields 



We start with free scalar field theory in two dimensions subject to periodic 
boundary conditions 



(j>{t,x - L) = (p{t,x + L). 



(1) 



In our discussion t is Euclidean time, analytically continued from imaginary 
Minkowskian time. Let J7 be an external source satisfying the same boundary 
conditions and which vanishes as \t\ —^ cx3. The Euclidean generating functional 
in the presence of J" is given by 

Z[J] ex Jv^c^p {-/_°°^dt/_\dx [i ((f )2 + (§1)2) + i^02 ^j^-^Y (2) 
We now expand in terms of periodic-box modes, 

These are also eigenmodes of momentum and each (^„ or »/„ carries momentum 
^. In terms of these modes, we have 



Z[J] « [UnT^^ne^p{-JZ,dtEr 
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For notational purposes we will define 






Using the Feynman-Kac formula, we find 



where 
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(4) 
(5) 
(6) 
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and |0) is the ground state of Ha. Since Hq is the usual equal time Hamiltonian, 
|0) is the vacuum. Ho consists of a set of decoupled harmonic oscillators, and 
it is straightforward to calculate the two-point correlation functions. 



(0|<^_„(t2)<^„(tl)|0) = 



SJ„{t2) SJ-^{ti) 



Z\J] 



J=0 



= ^o^v[-u:^\t2-t^\. (8) 



We now include a (j)'^ interaction term as well as a counterterm Hamiltonian, 
which we denote as Hc.t.- The new Hamiltonian is 



Hj=En 
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We will regulate the sums over momentum modes by choosing some large pos- 
itive number A^max and throwing out all high-momentum modes (/„ such that 
|n| > -/Vmax- This corresponds to a momentum cutoff 



A2 



' N„ 



(10) 



In two-dimensional (/> theory, renormalization can be implemented by normal 
ordering the 0^ interaction term. This corresponds to cancelling diagrams of 
the type shown in Figure 1. Using (g|), we find 
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where 

For the remainder of our discussion we will use the Hamiltonian 
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3 The 02 phase transition 



The existence of a second-order phase transition in two-dimensional 4>'^ theory 
has been derived in the literature Q [Q . The phase transition is due to devel- 
oping a non-zero expectation value and the resultant spontaneous breaking of 
— > — reflection symmetry. It is generally believed that this theory belongs 
to the same universality class as the two-dimensional Ising model and there- 
fore shares the same critical exponents. In this section we apply periodic field 
methods to (j>'^ theory in order to determine the critical coupling and critical 
exponents v and (3. We have chosen v and (3 since these are, in our opinion, 
the easiest to determine from direct computations. All other exponents can 
be derived from these using well-known scaling laws. From the Ising model 
predictions we expect 



^=1, /? = 



(14) 



v is the exponent associated with the inverse correlation length or, equiva- 
lently, the mass of the one-particle state. We will determine the behavior of the 
mass as we approach the critical point from the symmetric phase of the theory. 



Let I a) be any state even under reflection symmetry. We consider the matrix 
element 

f{t) = {a\qoexpi-tH)qo\a). (15) 

Inserting energy eigenstates \i) satisfying H \i) = Ei \i), we have 

/W = E.expM£;.OI(*ko|a)|'. (16) 

Since \a) and |0) are even under reflection symmetry and go is odd, the vacuum 



contribution to the sum in (16) vanishes. In the Hmit i ^ oo, (16) is dominated 
by the contribution of the next lowest energy state, the one-particle state at 
rest.^ In this limit we have 

f{t) ^ e-™*, (17) 

where m is the mass. We can compute f{t) numerically using the method of 
diffusion Monte Carlo (DMC) . The idea of DMC is to model the dynamics of the 
imaginary-time Schrodinger equation using the diffusion and decay/production 
of simulated particles. The kinetic energy term in the Hamiltonian determines 
the diffusion rate of the particles (usually called replicas) and the potential 
energy term determines the local decay /production rate. A self-contained in- 
troduction to DMC can be found in j|] . 

V is defined by the behavior of m near the critical coupling Ac, 

m^{K~XY. (18) 

Once we determine f{t) using DMC simulations we can extract m and v using 
curve-fitting techniques. We have calculated m as a function of A for several 
different values of L and A'max- Results from these calculations are presented 
in the next section. 

(3 is the critical exponent describing the behavior of the vacuum expectation 
value. In the symmetric phase the vacuum state is unique and invariant under 
the reflection transformation cj) ^ —(p (or equivalently q„ -^ —qn, for each n). 
In the broken-symmetry phase the vacuum is degenerate as L — > oo, and qo, the 
zero-momentum mode, develops a vacuum expectation value. In the L —> oo 
limit tunnelling between vacuum states is forbidden. One ground state, |0+) , is 
non-zero only for values qo > and the other, |0~), is non-zero only for qo < 0. 
We will choose |0+) and |0~) to be unit normalized. 

Let us now define the symmetric and antisymmetric combinations, 

I0') = 72(|0^) + |0")) (19) 



•^We are assuming that this contribution does not also vanish. This is generally true, and 
we can always vary \a) to make it so. 



We will select the relative phases of |0~) and |0"'") so that they transform from 
one to the other under reflection symmetry. [O**) and |0°) are then symmetric 
and antisymmetric (respectively) under reflection symmetry. 

To avoid notational confusion in the following, we will write g„ to denote the 
quantum-mechanical position operator and q„ for the corresponding ordinary 

variable. Let IVz} be the spectral family associated with the operator 

^ J z^ ( — 00,00) 

-^.n This implies that / dVz is a projection operator whose action on a general 
wavefunction ^ is 

(^fjVz)^{qo,qi,---)^0{^-a)9{b^^)^{qo,q^,---). (20) 

From the support properties of |0^) and |0~), we deduce 

|0+) - V2j^dVz |0^) = V2j^dVz lO"^) (21) 

Using our new spectral language, we can write 

-4^^ r zdVz. (22) 

,/2L J-00 '2: \ / 

We now consider the vacuum expectation value (0+| 4> |0+).q] Making use of 
translational invariance, we have 

(0+1 d^ |0+) = ^S-Ldx (0+1 </-(t,:r) |0+) - (0+1 Jj |0+) . (23) 

From ( pi| ) and (£2), we conclude that 

{Q+\<l>\Q+)=2mj^zdVzm^2J^dzzg{z), (24) 

where 

9{z) = m^m. (25) 

g{z) satisfies the normalization condition 

ST^dzg{z)^2S^dz9{z) = l. (26) 

In our calculations wc will be working with large but finite L. In this 
case the ground state degeneracy is not exact, and the symmetric state |0'') is 
slightly lower in energy than the antisymmetric state |0°). We can now use 
this observation (that |0*) is the lowest energy state) to rewrite g[z) as 

a(7\ - /n"! '^^^ W\ - Utt, f'\^^'P{-t"}^^^v{-^H}\b) 

9\^) - \" I "dT 1^ / - /i,™ (b|Toxp{-2tff}|f,) ' 1^0 



*The extra factor . — has been included for later convenience. 

^We could also consider (0~ I </< |0~). By reflection symmetry (0~| |0~) = — (O"*"! </> 0"*") . 



where \b) is any state such that (0^|6) is non-zero. 
For free field theory g(z) can be exactly calculated, 



g[z) = ^^e"^^^^ . (28) 



For non-trivial coupling we can calculate the right-hand side of 127^) using DMC 
methods. In Figure 2 we have plotted g{z) for L = 2.5tt and L = Stt. In 
each case |j = 2.76 and A = 4. All quantities are measured in units where 
^ — 1. As can be seen, the distributions are bimodal and the maxima for both 
curves occur near ±0.55. We observe that the peaks are taller and narrower 
for larger L. This is consistent with our intuitive picture of fluctuations in the 
functional integral. For a small but fixed deviation in the average value of 
(j), the net change in an extensive quantity such as the action or total energy 
scales proportionally with the size of system. Consequently the average size of 
the fluctuations must decrease with L. We can estimate the amplitude of the 
fluctuations, Acj), by assuming a quadratic dependence in A0 about the local 
minimum. The net effect of the fluctuation should not scale with i, and we 
conclude thatQ 

/ \ ^spatial dim. 

"^^-{vl) -Vl- (29) 

This agrees with the free field result in (OS) and also appears to be consistent 
with the peak widths plotted in Figure 2. 

Let Zmax be the location of the non-negative maximum of g{z). Since g{z) 
becomes sharply peaked as L -^ oo, 

(O+l (j) |0+) = 2j^dz zg{z) « 2z,„ax/o°°rf2g(z) = z^ax- (30) 

This gives us another option for calculating the vacuum expectation value. We 
can either integrate 2zg{z) or read the location of the maximal point z^ax- 
Both will converge to the same value as L — > oo. However, the Zmax result is 
less prone to systematic error generated by the 0{-j=) fluctuations described 

above. [| We will therefore use 

(0+|(/)|0+)-Zmax. (31) 

P is defined by the behavior of the vacuum expectation value as we approach 
the critical coupling, 

(0+|(/.|0+)^(A-A,)'3. (32) 

Using DMC methods, we have computed the A dependence of the vacuum ex- 
pectation value for several values of L and A^max- The results are shown in the 
next section. 



^The dimension of time does not enter here since we are considering properties of the 
vacuum, the ground state of the Hamiltonian defined at a given time. This is in contrast 
with lattice calculations which usually consider the quantity (<^), the average of <f> over all 
space and time. 

'^We can see this explicitly in free field theory, where the vacuum expectation value should 
vanish. Zmax = as desired, but 2jr°dz zg{z) = . 



4 Results 

The results of our diffusion Monte Carlo simulations are presented here. For 
each set of parameters L and iVmax, the curves for m and (0+| |0^) near the 
critical coupling have been fitted using the parameterized forms 



«(¥-^)' (33) 



and 

(0l^|0+)=5(A-4;)^ (34) 

For the mass curves, data points with m < L~^ have correlation lengths ex- 
ceeding the size of our system and are of questionable significance. We have 
therefore fit these curves two different ways, once using all data points and a 
second time excluding small m values. The curves for the data set L = 2.57r 
and A^max = 10 are shown in Figures 3 and 4.^ The error bars represent an 
estimate of the error due to Monte Carlo statistical fluctuations and, in the case 
of the mass data, pollution due to higher energy states. Let Xd denote the 
reduced chi-squared value for d degrees of freedom. The results of the fits are 
as follows: 
L = 2.57r, Mnax = 8 (A = 3.2): 

^ = 3.1 ± 0.2, J/ = 1.1 ± 0.1, a = 0.31 ± 0.02, x?5 = 0.81 (all), (35) 

^ = 3.1 ± 0.2, 1/ = 1.0 ± 0.1, a = 0.30 ± 0.03, Xs = 0.96 (partial), 

^ = 2.5 ± 0.1, 13 = 0.15 ± 0.02, b = 0.65 ± 0.05, xl = 0-45 

L = 2.57r, 7V„,ax = 10, (A = 4): 

^ = 2.9 ± 0.2, 1/ = 1.1 ± 0.1, a = 0.35 ± 0.02, x?2 = 0-46 (all), (36) 

^ = 2.9 ± 0.2, 1/ = 1.0 ± 0.1, a = 0.37 ± 0.03, xl = 0.37 (partial), 

4^ = 2.5 ± 0.1, [3 = 0.18 ± 0.01, b = 0.70 ± 0.03, Xg = 0.29 

L = 5^, 7V„,ax = 16, (A = 3.2): 

^ = 3.0 ± 0.2, 1/ = 1.2 ± 0.1, a = 0.32 ± 0.03, xlz = 0-87 (all), (37) 

^ = 3.1 ± 0.2, ;/ = 1.2 ± 0.1, a = 0.30 ± 0.03, x?o = 0-88 (partial), 

^ = 2.7 ± 0.1, [3 = 0.12 ± 0.02, b = 0.62 ± 0.04, xl = 0.50 

L = 5^, TVmax = 20, (A = 4): 

^ = 2.8 ± 0.2, t/ = 1.2 ± 0.1, a = 0.35 ± 0.06, x?i = 1-2 (all), (38) 

^ = 2.9 ± 0.2, v=\.i± 0.1, a = 0.36 ± 0.06, Xs = 1-1 (partial), 

^ = 2.5 ± 0.1, 13 = 0.11 ± 0.02, b = 0.65 ± 0.04, xl = 0.88. 



'As mentioned before, we are using units where ^ = 1. 



These results are subject to errors due to the finite size L and finite cutoff scale 
A. We will use our data for different values of L and A to extrapolate to the 
limit L — > oo, A -^ oo. For the parameters i/, a, /3, h we use the naive asymptotic 
form 

x{L, K) = x + ^.TA2 + ■^xl2 +■■■ . (39) 

For the critical couplings A™ and Ac , however, we modify the finite L correction 
according to the finite-size scaling hypothesis H 

A(L,A) = A+^Aa. + ^A|^|V. + •••• (40) 

The results we find aren 



^ = 2.5 ± 0.2 ± 0.1, 1/ = 1.3 ± 0.2 ± 0.1, a = 0.43 ± 0.05 ± 0.02 (all), (41) 
^ = 2.3 ± 0.2 ± 0.1, v = lA± 0.2 ± 0.1, a = 0.48 ± 0.06 ± 0.02 (partial), 
^ = 2.5 ±0.1 ±0.1, /3 = 0.13 ±0.02 ±0.01, 6 = 0.71 ± 0.04 ± 0.03. 

The first error bounds include inaccuracies due to Monte Carlo statistics, higher 
energy states (for the mass curves), and extrapolation. The second error bounds 
represent estimates of the systematic errors due to our choice of initial state, 
time step parameter, and bin sizes in the DMC simulations. For data generated 
from the mass curves, the main source of error was due to extrapolation. The 
extrapolation error for the vacuum expectation value data was still the most 
significant, though considerably smaller than that for the mass calculations. 
This reduction is probably a result of the method used to measure the vacuum 
expectation value.^ Our results for the critical exponents are consistent with 
the Ising model predictions 
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The results for the critical coupling 
recently obtained lattice result R]P] 



5 Summary 

We have discussed the generalization of spherical field theory to other modal 
expansion methods, in particular, periodic field theory. Using periodic field 

^Due to the relatively weak dependence on L and A, the reduced chi-squared values for 
our extrapolation fits are quite small (<S 1) and do not provide a useful statistical measure. 

-'"We are refering to the result ^0+ 1 </< 1 0+ ) = Zmax, which eliminates peak broadening effects 
for finite L. This is discussed at the end of the previous section. 

^^ Critical exponents were not measured in this study. 



methods we have analyzed two-dimensional ^^ theory and computed the critical 
coupling and critical exponents v and /3 associated with spontaneous breaking of 
— > — (/) reflection symmetry. Our value of the critical coupling is in agreement 
with a recent lattice calculation, and our values for the critical exponents are 
consistent with the critical exponents of the two-dimensional Ising model. This 
lends support to the popular belief that the two theories belong to the same 
universality class. 

The full set of diffusion Monte Carlo computations used in our analysis 
required about 30 hours on a 350 MHz PC processor. Complete codes can 
be obtained upon request from the authors. The required computational time 
appears to be dominated by the number of operations required to update the 
Hamiltonian, which scales as A^max- Errors can be reduced quite substantially 
by using larger values of L and iVmax and utilizing large-scale parallel processing. 
No less important, however, is that periodic field theory provides a simple and 
efficient approach to studying non-perturbative phenomena with only modest 
computer resources. Improvements are now under way to utilize fast Fourier 
transform methods ||] and increase the computational speed. Future studies 
have been planned to analyze phase transitions in other field theory models. 
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Figures 



Figure 1. The only divergent diagram, which can be cancelled by normal order- 
ing. 

Figure 2. Plot of g{z) for L = 2.5tt and L = Stt. In each case ^ ~ 2.76 anc 
A = 4. 



Figure 3. Plot of to as a function of ^ for i = 2.57r, A^max = 10. 
Figure 4. Plot of (0+| (f> |0+) as a function of ^ for i = 2.5tt, A^max = 10 
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Figure 1 
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Figure 3 
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